Suprathreshold stochastic resonance: suppression of noise by strong periodic signals 
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We consider an overdamped Brownian motion in "quartic" potential subjected to periodic driving. 
This system for the case of a weak periodic driving has been intensively studied during past decade 
within the context of stochastic resonance. We had demonstrated, that for the case of predominantly 
suprathreshold driving the noise in the output signal is strongly suppressed at certain frequency 
range: the signal-to-noise ratio demonstrates resonant behavior as function of frequency. 

The switching dynamics of an overdamped brownian particle in a "quartic" potential has been intensively studied 
in the past decade in the frame of stochastic resonance (SR) phenomenon Stochastic resonance is a nonlinear 
noise-mediated cooperative phenomenon wherein the coherent response to a deterministic signal can be enhanced in 
the presence of an optimal amount of noise. It had been observed in a wide variety of electronic systems, such as, 
e.g., lasers §, Schmitt triggers |], tunnel diodes § and SQUIDs @, §. It is known, that for a singlelevel threshold 
systems SR has been observed for the case of weak (underthreshold) driving (although residual SR effects are 
known to occur for marginally suprathreshold signals M). In this case as well as in the case of suprathreshold SR in 
multilevel threshold systems the manifestation of SR is resonant behavior of signal-to- noise ratio {SNR) or other 
relevant characteristics as function of noise intensity. 

However, most of practical devices are operating in the predominantly suprathreshold regime, when the transition 
from one state to another one over potential barrier occurs deterministically and noise is the only disturbing factor 
leading to erroneous switching. As an example of such a situation we can refer to microwave hysteretic SQUID 
1^^ , whose dynamics is described by the model of Brownian motion in bistable potential subjected to strong periodic 
driving at the given frequency. Such SQUID represents a clear example of the device, by its very basic idea operating 
in strongly suprathreshold regime (its important characteristic is a function of amplitude of output signal versus the 
amplitude of input driving) and for which presence of noise leads to earlier transitions over potential barrier that 
results in error of the measured dc magnetic flux. For microwave hysteretic SQUID long standing problem is known: 
at which parameters it should operate in order to demonstrate maximal sensitivity. On one hand it is known, that 
with increase of pumping frequency the sensitivity of the SQUID should improve, on the other hand it is known that 
at frequencies higher than the cut-off frequency, the performance of the device should degrade. Some qualitative 
treatment of this problem has been done in the past 0, [0, but the question is still opened. 

Few papers, where the case of strong periodic driving was considered both in classical [ p^ , and quantum 
systems p^ , were addressed to investigation of area of hysteretic loop and its resonant behavior as function of 
frequency of driving signal was demonstrated. But, to our knowledge, the investigation of signal-to-noise ratio as 
function of frequency was not performed for the case of strongly suprathreshold driving. This may be explained by 
the fact, that the most studies were restricted by adiabatic approximation, where frequency dependence of SNR could 
not be investigated in detail or by linear response theory, where the driving amplitude was supposed to be small. 
Nevertheless, in the frame of linear response theory some weak resonant frequency dependence of signal-to-noise ratio 
was recently observed for a particular case of piecewise rectangular potential . 

Recently we have shown using description via temporal characteristics, that in a dynamical system with noise 
subjected to a strong periodic driving, significant suppression of noise is possible in a certain frequency range p^ ]. 

The aim of the present paper is to consider fluctuational dynamics of a Brownian particle in a " quartic" potential, 
study the signal-to-noise ratio and demonstrate that for the case of a strong (suprathreshold) periodic driving the 
SNR has resonant behavior as function of frequency. 

Consider a process of Brownian diffusion in a potential profile 



U{x,t) = te'* - ax'^ -t- xAsm{nt + ip), 



(1) 



where Lp is initial phase. It is known that the probability density W{x, t) of the Brownian particle in the overdamped 
limit (Markov process) satisfies the Fokkcr-Planck equation (FPE) 
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Here G(x,t) is the probability current, B = h/kT, h is the viscosity (in computer simulations we put h = 1), T is the 
temperature, k is the Boltzmann constant and u{x,t) = U{x,t)/kT is the dimensionless potential profile. The initial 
and the boundary conditions have the following form: 
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W{x, 0) = 5{x - xq), G(±c5o, t) = 0. 



(3) 



In computer simulations we had chosen the following parameters of the potential: b = 1, a = 2. With such a choice 
the coordinates of minima equal Xmin = il, the barrier height AU = 1, the critical amplitude Ac is around 1.5 
and we have chosen A = 2 to be far enough from Ac- We note, that this is not the case "just above the threshold 
level", considered in 0, but indeed strong driving: the amplitude A — 2 is far above the dynamic threshold. We also 
performed the analysis for ^ = 3; 4; 5: the results are qualitatively the same, only SNR rises accordingly. 

The quantity of our interest is the SNR. In accordance with we denote SNR as: 

SNR = ^ \im^ I S{uj)duj, (4) 

where 



O-A, 



+ OC 



S{uj) = j e-"^^K[T]dT (5) 

— OO 

is the spectral density, Sm{S^) is noisy pedestal at the driving frequency and K[t\ is the correlation function: 

K[T] = {{x{t + T)xm. (6) 

where the inner brackets denote the ensemble average and outer brackets indicate the average over initial phase ip. 

In order to obtain the correlation function K[t\ we solved the Eq. (|^) numerically, using the Crank-Nicholson 
scheme. 

In Fig. 1 the spectral density S{uj) is presented for kT = 0.1 (the delta-spikes at 1-st and some higher harmonics 
are out of figure in order to enhance the noise part). One can see that the form of S{lo) significantly depends on 
driving frequency fi, while the amplitude of the output signal is monotonically decreasing function of fi. In order 
to study the resonant behavior of spectral density, let us plot the SNR as function of driving frequency f2. From 
Fig. 2 one can see, that SNR as function of has strongly pronounced maximum. The location of this maximum 
at = ^max approximately corresponds to the time scale matching condition: 

max ~ "^/^min^ whcrC T^^^ is the 

minimal transition time from one state to another one. The existence of optimal driving frequency may be explained 
in the following way. Let us consider the case of adiabatically slow driving. If noise is absent, the escape would 
occur only after the corresponding potential barrier would disappear. If we add some small amount of noise, the 
escape would occur earlier than in the deterministic case at some nonzero barrier height. If the driving frequency 
is increased, the potential barrier height will decrease faster and the escape will occur at lower barrier, that is more 
close to the deterministic case. In the case, when the driving frequency is higher than the cut-off frequency of the 
system, D, > flc (where flc has a dynamical sense: for a given amplitude of the signal the performance of the system 
significantly degrades above certain frequency), the particle would never escape over potential barrier in the absence 
of noise and will remain in the vicinity of initial potential minimum, since there is not enough time to reach the basin 
of attraction of another state. Therefore, there is some frequency range, where at the given small noise intensity the 
escape will occur over the smallest potential barrier and namely in this case noise has minimal effect on the system 
that results to the maximal SNR. 

In the case when the driving frequency is higher than the cut-off frequency, > adding some amount of noise 
will help the particle to move to another state and the conventional stochastic resonance may be observed (see the 
inset of Fig. 2 for il = 1). Let us note, that we were not able to observe such intuitively obvious thing as increase 
of SNR for kT ^ in the case il > ftc (in this case for kT = we will have small but nonzero oscillations near 
the initial potential minimum) due to divergence of the algorithm for kT < 0.02 (fi = 1) and the last three points, 
presented in the inset of Fig. 2, correspond to kT — 0.1; 0.05; 0.02, respectively. 

In conclusion we have shown that in the dynamical system with noise driven by a strong sinusoidal signal (pre- 
dominantly suprathreshold driving) the influence of noise is significantly reduced in a certain frequency range: the 
signal-to-noise ratio is a resonant function of frequency of the driving signal. This effect is of real importance for 
applications since it may allow to operate a concrete device (like a SQUID or other electronic devices) in the regime 
of minimal noise-induced error. 

For driving frequencies higher than the cut-off frequency the conventional (but suprathreshold) stochastic resonance 
has been observed. 
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FIG. 1. Spectral density S{uj) with enhanced noise part, kT = 0.1, A = 2. 
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